NASA/TM— 2009-216096 



Analysis of Electromagnetic Wave Propagation in a 
Magnetized Re-Entry Plasma Sheath Via the 
Kinetic Equation 


Robert M. Manning 

Glenn Research Center, Cleveland, Ohio 


December 2009 


NASA STI Program ... in Profile 


Since its founding, NASA has been dedicated to the 
advancement of aeronautics and space science. The 
NASA Scientific and Technical Information (STI) 
program plays a key part in helping NASA maintain 
this important role. 

The NASA STI Program operates under the auspices 
of the Agency Chief Information Officer. It collects, 
organizes, provides for archiving, and disseminates 
NASA’s STI. The NASA STI program provides access 
to the NASA Aeronautics and Space Database and 
its public interface, the NASA Technical Reports 
Server, thus providing one of the largest collections 
of aeronautical and space science STI in the world. 
Results are published in both non-NASA channels 
and by NASA in the NASA STI Report Series, which 
includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or theoretical 
analysis. Includes compilations of significant 
scientific and technical data and information 
deemed to be of continuing reference value. 
NASA counterpart of peer-reviewed formal 
professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies that 
contain minimal annotation. Does not contain 
extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, often 
concerned with subjects having substantial 
public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific and 
technical material pertinent to NASA’s mission. 

Specialized services also include creating custom 

thesauri, building customized databases, organizing 

and publishing research results. 

For more information about the NASA STI 

program, see the following: 

• Access the NASA STI program home page at 
http://www. sti. nasa.gov 

• E-mail your question via the Internet to help@ 
sti.nasa.gov 

• Fax your question to the NASA STI Help Desk 
at 443-757-5803 

• Telephone the NASA STI Help Desk at 
443-757-5802 

• Write to: 

NASA Center for AeroSpace Information (CASI) 
7115 Standard Drive 
Hanover, MD 21076-1320 


NASA/TM— 2009-216096 



Analysis of Electromagnetic Wave Propagation in a 
Magnetized Re-Entry Plasma Sheath Via the 
Kinetic Equation 


Robert M. Manning 

Glenn Research Center, Cleveland, Ohio 


National Aeronautics and 
Space Administration 


Glenn Research Center 
Cleveland, Ohio 44135 


December 2009 


This report contains preliminary findings, 
subject to revision as analysis proceeds. 


Level of Review : This material has been technically reviewed by technical management. 


Available from 


NASA Center for Aerospace Information 
7115 Standard Drive 
Hanover, MD 21076-1320 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 


Available electronically at http://gltrs.grc.nasa.gov 


Analysis of Electromagnetic Wave Propagation in a Magnetized 
Re-Entry Plasma Sheath Via the Kinetic Equation 

Robert M. Manning 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

Based on a theoretical model of the propagation of electromagnetic waves through a hypersonically 
induced plasma, it has been demonstrated that the classical radioffequency communications blackout that 
is experienced during atmospheric reentry can be mitigated through the appropriate control of an external 
magnetic field of nominal magnitude. The model is based on the kinetic equation treatment of Vlasov and 
involves an analytical solution for the electric and magnetic fields within the plasma allowing for a 
description of the attendant transmission, reflection and absorption coefficients. The ability to transmit 
through the magnetized plasma is due to the ‘magnetic windows’ that are created within the plasma via 
the well-known ‘whistler modes’ of propagation. The case of 2 GHz transmission through a re-entry 
plasma is considered. The coefficients are found to be highly sensitive to the prevailing electron density 
and will thus require a dynamic control mechanism to vary the magnetic field as the plasma evolves 
through the re-entry phase. 

1.0 Introduction 

Electromagnetic wave (EM) propagation through a flowing plasma layer to maintain communications 
and navigation to hypersonic and space re-entry vehicles has been a problem for over 45 years. The 
classical radio frequency (RF) blackout results when the associated EM wave is reflected and absorbed by 
the free electrons that make-up the plasma sheath which envelops the re-entering vehicle. The plasma is 
the result of extreme heating of air by the strong shock wave that is created by the leading edges of the 
vehicle [Bletzinger, et al., 2005 and references therein]. Such a plasma has associated with it a 
characteristic frequency, the ‘plasma frequency’, below which RF transmission is reflected and absorbed; 
RF radiation above this frequency will easily pass through the plasma. Since this plasma frequency is 
proportional to the electron concentration within the plasma, this frequency will vary according to the 
altitude and shape of the vehicle, as well as the angle of attack. Thus, for a given communications 
scenario, a particular frequency is assigned; As the vehicle begins its re-entry into the earth’s atmosphere, 
shock waves are created and a plasma begins to be formed. At a particular electron concentration, the 
plasma frequency will exceed that of the communications link and transmission to and from the vehicle 
will cease. As the vehicle begins to decelerate, the electron density and thus the plasma frequency 
decrease and RF communication to and from the vehicle once again becomes possible. 

Several approaches have been advanced in order to mitigate RF blackout. The obvious ones are to use 
communications frequencies higher than that of the plasma frequency that will develop during re-entry. 
However, as mentioned above, a specific maximum plasma frequency is dependent upon several 
parameters and becomes very difficult to establish. A high power can be used to overcome the reflection 
and absorption of the plasma. This, however, severely complicates the design and operation of a 
communications system, especially the equipment on the re-entering vehicle. Another method is to inject 
electrophilic substances (quenchants) into the plasma flow field of the sheath so as to de-ionize the 
plasma (at least, severely reduce the electron concentration) and lower the plasma frequency below that of 
the communications link. An application of this method occurred during the re-entry of Gemini 3 in 1965 
[Schroeder and Russo, 1968] where it was demonstrated that plasma reflection and absorption can be 
significantly reduced to allow for re-entry communications. Another technique that can be employed is 
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that of establishing ‘magnetic windows’ within the plasma. Here, a static magnetic field is applied to the 
plasma to essentially establish ‘whistler modes’ of propagation [Usui, et al., 2000]. Early experiments 
using this concept were performed in 1964 [Russo and Hughes, 1964]. Here, a magnetic field of 750 
Gauss was used in a ground experiment employing a plasma from a solid rocket motor. It was shown that 
a signal improvement of about 20 dB can be realized. However, it was recently argued [Sharkey, 2003] 
that magnetic fields on the order if 10 4 Gauss are needed to penetrate a re-entry plasma. 

The ‘magnetic windows’ concept is very attractive in that, unlike the use of electrophilic substances 
where a supply of such material is required throughout the blackout period (which, for some planetary re- 
entry scenarios, can have a duration of 10 to 20 minutes), the magnitude of the static magnetic field can 
be adjusted as conditions require to maintain a window in the otherwise RF opaque plasma. The 
electrophilic technique is an ‘active’ mitigation method whereas the magnetic windows technique is a 
‘passive’ method. If it can be demonstrated that magnetic fields of nominal strength can be used to elicit a 
magnetic window within a plasma with re-entry parameters, and, additionally, it can be shown how the 
magnitude of this magnetic field needs to be adjusted to maintain a window as the plasma properties 
evolve, then the use of an external magnetic field as a mitigation technique will be a viable one. 

This prescription can only be accomplished by realistically modeling the plasma propagation 
environment. It is the purpose of this report to analytically study such a situation. A mathematical model 
will is constructed, based on the Vlasov equations, i.e., the Maxwell equations supplemented with the 
kinetic (Boltzmann) equation describing the electron distribution, that will capture the propagation 
process through a flowing plasma immersed in an external magnetic field. The model will then be solved 
in the approximation of weak spatial dispersion; the case of strong spatial dispersion will be treated in a 
forthcoming report. It is shown that an applied magnetic field of nominal strength can alter the plasma 
and create “magnetic windows” through which electromagnetic radiation can propagate. 

2.0 The Initial Equations 

2.1 Incorporating the Maxwell Equations with Kinetic Theory — The Vlasov Equations 

The basic starting point for modeling electromagnetic wave propagation through a flowing 
hypersonic plasma are, of course, the Maxwell Equations 


V ■ E = 4jrp 

(2.1) 

V- 5 = 0 

(2.2) 

VxE + i— = 0 

(2.3) 

c 8t 

- - 1 dE 4n - 

(2.4) 

V x 6 — — = —j 
c 8t c 


The expressions relating the charge density p and current density j within the plasma are given in terms 
of the statistical distribution functions f s (r,v,t) governing the charge carriers of the species s at 
position r and velocity 6 at time t 


oo 

P = ^e* n s { fs (r,u,t)d 3 v 

s —oo 


(2.5) 
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( 2 . 6 ) 


j = ^ e s n s { uf s (r,u,t)d 3 v 

S -00 

Here, 5 denotes either electron or ion charge carriers where n s = n ion is the number of positively charged 
ions, e s = e ion = +e or n s = n e is the number of negatively charged electrons, e s = e e =-e. The statistical 
distribution function for species s f s (r,u,t) is given by the Boltzmann equation incorporating collisions 
via the Krook model [Tanenbaum, 1967; Bhatnagar, et al., 1954] for collisions, viz, 

^-+6 - Vrfs + Fs ■ Vo/ = -v, (/, - fs 0 ) (2.7) 

where v s is the effective collision frequency and f s o{r,u) is the initial equilibrium distribution of the s- 
th species. The model for the effective collision frequency used here, along with the basis of the Krook 
model, is discussed in Appendix A. The force entering Equation (2.7) is given by the Lorentz force 
involving both the electric E and magnetic B fields of the wave as well as an externally applied 
magnetic field 
B 0 , 

F s =^(e + -x(b + B 0 ) ) (2.8) 

m s \ c ) 

Collectively, the Maxwell equations, Equations (2.1) to (2.4), the Equations (2.5) and (2.6), and the 
Boltzmann equation, Equation (2.7) with Equation (2.8) are known as the Vlasov equations. 

These coupled integral and differential equations must now be applied to the situation of a 
hypersonically flowing plasma upon which is incident a plane electromagnetic wave. The situation is 
depicted in Figure 1 . The tangential velocity V T is due to the plasma flowing along the hypersonic vehicle 

surface. The velocity V L is due to the motion of the hypersonic vehicle toward the observer. The plasma 

is taken to be of infinite extent along the x and y axes. An incident plane electromagnetic field impinges 
on the moving plasma layer from the transmitter fixed on the earth and travels along the z-axis in the -z 
direction. The constant homogeneous magnetic field emanates from an antenna on the surface of the 
vehicle along the z-axis in the +z direction. The surface in the jc -y plane is taken to be along the solid 
surface of the vehicle and that at the front of the layer of thickness L is taken to be open to the 
atmosphere. The values of the E and B fields of the electromagnetic wave are reckoned with respect to 
the observer in the rest frame with respect to the moving plasma layer. Given these conditions, one can 
write for the functional dependence of the field and current density, 

E = E(z,t)=E x (z,t)x + E y (z,t)y (2.9a) 

B = B(z,t)= B x {z,t)x + B y {z,t)y + Bqz (2.9b) 

j = j{ z ,t)= j x {z, t)x + j y (z, t)y (2.9c) 

The problem is most easily dealt with by transforming the fields into the reference frame moving with the 
plasma. 
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2.2 Transformation of the Fields Into the Reference Frame Moving with the Plasma 

The transformation of the fields E and B in the reference frame of the observer (transmitter) to 
those El and B' that are seen in the frame moving with the plasma layer is given by the well-known 
relations [Jackson, 1975] 


II 

+ 

■C30[ 

X 

1 

+ 1, 

(2.10) 

El = y(s' Px£) Y p(j3-s) 

(2.11) 

where 



(2.12) 

and c is the velocity of light. Also, for the current density 

/ = y(i-Pp) 

(2.13) 

The invariance of the Maxwell equations has that, in the frame moving with the plasma, 


V'x£' + -— = 0 
c df 

(2.14) 

1 dE' 4 k 

V'xS 1 = — f 

c di c 

(2.15) 


where V' and t' are the appropriately transformed divergence and time. Consider now the extreme case 
in which a hypersonic velocity of Mach 35 is realized, i.e., Vj = Vi « 1.2x 10 4 m/s. This yields 
P ~ 4.0x 10 -5 « 1 . In this case, the transformed equations Equations (2.13) and (2.14) essentially reduce 
to those in the rest frame, i.e., Equations (2.3) and (2.4). Similarly, for the homogeneous applied magnetic 
field, B'q ~ Bq. Hence, to within a first order approximation, Equations (2.1) to (2.4) can be employed in 
the reference frame moving with the plasma. One can now incorporate these relations in the evaluation of 
Equation (2.7). 

One should also address the Doppler shifts that will be incurred in the transformation to the reference 
frame moving with the plasma. This will not be addressed at this point as it is the goal of this work to 
establish the possibility of the magnetic windows concept. A more careful study with will appear in a 
forthcoming work will incorporate Doppler effects. 

The final goal in the analysis that follows is the evaluation of the reflection, transmission, and 
absorption coefficients of the plane wave interacting with the plasma layer; the ‘flow chart’ of this 
process is depicted in Figure 2. The scattering process, as well as its solution, goes as follows. The 
electric and magnetic fields of the plane electromagnetic wave induces a Lorentz force on the plasma 
electrons. The induced force, in turn, produces a variation in the electron distribution function governing 
the position and momentum of the electrons. This gives rise to an induced current within the plasma that 
produces an associated variation of the fields within the plasma which, once again, induces an additional 
Lorentz force. This non-linear process finally results in the establishment of field distributions at the 
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boundaries of the plasma that determines the corresponding reflection, transmission, and absorption 
coefficients of the incident field. 


3.0 Application of the Foregoing to the Boltzmann Equation 

3.1 Reduction of the Equation 

Since the current density is only a function of the longitudinal coordinate z, so too will be the 
distribution function of Equation (2.7) 

/s(^C,f)=/ s (z,C,f) (3.1) 


The associated equilibrium distribution function is taken to be given by the Maxwell distribution 


fsO («) = fso (|u|) = 

\n j 


\i /2 r i s 3 


V Ts ) 


exp 


'77T 

V V Ts J 


(3.2) 


where the thermal velocity of the s-th species is given by 


UTs = 


UbTs 

m s 


(3.3) 


Here, k B is Boltzmann’s constant, T s is the absolute temperature of the 5 -th species of charges with mass 
m s . The distribution function f s (z,v,t ) will be taken to be related to that of Equation (3.2) by a small 
perturbation ^(zjU,*), viz, 

/j(z,6,f)= f s0 jo|)+c|)(z,o,t),c|)(z,o,t)« / s0 (jo|) (3.4) 

In order for (f), 5 (z,o,i) to be treated as a perturbation, the fields of the incident wave must also be treated 
as a first order perturbation; writing the temporal dependence of the fields of the incident electromagnetic 
wave as 


E x (z,t) = E x (z)exp(-i(ot), etc. 


one has from Equations (2.9a) and (2.9b) 


E = (E , x (z)jc + E' >) (z))exp(-ia)f) 

(3.5) 

B = [B x (z)x + B y (z)yjexp(-z'fi)t) + Bq (z) 

(3.6) 


in which one must now require for a perturbation solution, E x ,E y ,B x ,B y « Bq . Finally, since the time- 

harmonic fields of Equations (3.5) and (3.6) are taken to be the source of perturbation of the distribution 
function given by Equation (3.4), one can further write 

§ s (z,o,f) = ^ (z,o)exp(-fcot) (3.7) 
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Substituting Equations (3.4) to (3.7) into Equation (2.7) and dropping quantities that are second order 
perturbations with respect to the fields and to ^(zjC) yields 


-/(co + z'Vj)^ (z,o)exp(-za>t) + o z ’ ^ exp(-z'g)t) 


+ 


'e, v 


\ m s J 


Ve/,o(|6|)+ 7 X *(*)■ Vo /,0 (|o|)jexp(-ioi)f) 

7 * B 0 z- V»f s0 (|o|)+— x B 0 z- (z,o)exp(-z'oot)|j = 0 
n C J 


+ 


m 


-s JL 


(3.8) 


However, remembering the identity (ox b)- d = 0 , one has that 

7 * 5(z)- Vufso (| 6 |) = ^x B(z} U 8fs0 


o So 


(3.9) 


and similarly for the term (o/c)x Bqz- V 5 / s q(|o|) . (The term (o/c)x B 0 z- (z,o) ^ 0 since c() i 5 (z,o) 
is not a function of the scalar V .) Hence, Equation (3.8) reduces to 


, \ / \ 5^(z,o) , , 

-z'(co + iv s jf, (z,o) + o z f ~exp(-za)t) 


Sz 


+ 


m 


( „ V 


Z(z)-V 0 ^(z,o)+ —x Bqz •V 0 ^(z,o) = 0 


s J 


V m s A ' 


(3.10) 


Defining the cyclotron frequency 




_ e^Bo 

m s c 


(3.11) 


and expanding ox z- Vo([), s . = (o ^ 9/9o x - o x djdv> y )t>, as well as E- V 6 f s0 = (E x d/dv x + E y d/dv y )f s0 , 
Equation (3.10) becomes 


o 2 ^-z(co-hv,)^ + o) c 


o. 


So, 


-o. 


So 


+ 


y J 


V 


\ m s J 


So, 


•+ E v 


So 


fsO = 0 (3-12) 


y J 


This differential equation is most easily dealt with in plane-polar coordinates Ox and 0 defined by 
o x = ox cos0,o^ = Ox sin0 . In this instance, one has 


3 a d 
- COS0- 


So, 


sin 0 S S . S cos 0 S 
= sin 0 + - 


Sox ox 90 So, 


Sox ox 90 


(3.13) 


allowing Equation (3.12) to be written as (noting that f s o is isotropic and independent of the polar angle 

9) 


9ct> 


Sz 


- - ® + i‘v s )<t>j. - (D Ci — + 


90 


2- |(^cos0 + ^sine)^. = O (3.14) 

j Sox 
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Finally, using the auxiliary fields [Gross, 1951; Bell and Buneman, 1964] defined by 

E l = E x — iEy ,E R = E x + iE y 


(3.15) 


E x cos0 + E y sin 0 = ~{El exp(i0)+ E R exp(- z'0)) 
Substituting this into Equation (3.14) and rearranging terms yields 


(3.16) 


tOca ^ ~ ^ ^ ° ~ + i ( co + iv s ) ( z,i o ) = (£ L exp(z0) + £* exp(-z0))|^. = 0 (3.17) 

(A 7 OZ ^ lWl s J 

This form of the differential equation can now be solved using the method of characteristics. Writing the 
equation as a relation along its characteristic curve, one has 

dti> s 5cL 50 56* dz 56. 56* 

= 1 = C0 r? V- 

dt 50 dt dz dt 50 dz 


Therefore, 


0 = (o cs t + 0',z = -v z t 


(3.18) 


which combines to give 0 = -(co CiS /u z )z + 0' . Changing back to the variables 0 andz, d/dt = -v z (d/dz) 
and Equation (3.17) becomes the first order differential equation 


h - + »(to + iv s (z,o) 


t E l exp(-z(ro Ci /o z )z + 0 ) + E R exp(z(ro CiS /u Z )z- z0' ))|^- = 0 


This equation must now be solved across the plasma layer of thickness L, 0 <z<L. 

3.2 Solution of the Differential Equation Across the Plasma Layer 

Consider the Fourier series of the function c|) s (z,u) over the interval 0 < z < L ; 


(3.19) 


<ta(z,6)= X ^( K /> u ) ex P(*‘ K / z )’ k /=T 

/=-o 0 L 


(3.20) 


If 

( K l O ) = — J ( z,\ u ) exp ( -IK/Z ). dz 


(3.21) 


Thus, multiplying Equation (3.19) by exp (-z/c ; z) and integrating over z, using Equation (3.21) gives 
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(3.22) 




l K/ 


iCco + iv.)! / . , ./ 

— -I L = ♦, (O,o)-^ {L , - 

j 


o- A 2m 




’S J 


} dfsO 

9ox 


El' 


K /+' 


CO. 


+ £*' 


K/-- 


co. 




J z yj 


where exp(-z'K/Z) = exp(-/7t) = (-1| and i?/ (V) = (/r)exp(z‘0), Ejf (/r) = E R (/r)exp(-z‘0) . At 

this point, contact must be made to the current density expression given by Equation (2.6); converting to 
plane polar coordinates, 


oo 2n oo 

l{ z ) = Yj e ° n s \ | |o x (xcos0+j)sin0)/ s (z,oj_,o z ,t)o 1 rfoj L rferfo z (3.23) 

s -oo 0 0 

As before, considering the identity 


xcos0 + ysin0 = -j[(jc-zy)exp(z0) + (Jc + zy)exp(-z0)J 


and using it in Equation (3.23) and taking the scalar product with the particular vector x + iy yields, 

oo 2n oo 

j R {z,t) = j x (z,t)+ijy{z,t) = Y 4 e s n s J j |v)iexp(z0)/ 5 (z,o ± ,o z ,t)c?o_ L c/0c/o z (3.24) 

S -00 0 0 

Thus, the current density has been converted to a complex scalar current density in the transverse plane. 
Substituting now Equations (3.4) and (3.7) into Equation (3.24), and using Equation (3.2), one as that the 
/so(|o|) term vanishes upon performing the u z integration. Hence, what survives is the result 


oo 2n oo 

jR{ z ) = 'Yj e s n s \ | jo^exp(z0)c|) i (z,o 1 ,o z )^o 1 c?0Jo z (3.25) 

s -OO 0 0 

where j R (z,t)= j R [z)ex p(- mt ) . Thus, the current density associated with the polarized wave fields E R 
and E l of Equation (3.15) are also polarized; however E R and E L remain uncoupled. Therefore, one can 
use either E R or E L separately in Equation (3.22). Selecting E R , i.e., the whistler mode E R = E X + iE y , 

Equation (3.22) becomes, upon using the definition kJ = K/ -co C5 /o z , 


where 




) = [(k',-Q 5 )zX] JA # (o.l,o. 



Q 


S 


CO-C Ocy+fVj 
^z 


(3.26) 


(3.27) 


and 
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A 5(ui,u z )=c|) i (0,u 1 ,u z )-c|) J (L,uj.,u z X-iy (3.28) 

Now, rewriting Equation (3.25) by limiting the u z integration in velocity space to the interval 0 < u z < oo 
gives 


00 231 00 

■M z ) = X e » a, 'f J Jwiexp(i0)O 5 (z,i) ± ,o z )c?o 1 rf0</o z (3.29) 

•S 0 0 0 

where 

(^,o ± ,o z ) 3 <t» s (z,o 1 ,o z ) + c|) J (z,o ± ,-o z ) (3.30) 

The Fourier transform of Equation (3.29) is, using the form of Equation (3.21), 

oo 2 tc oo 

h{ K i) = ^ e s n s\ I J exp(i6)d> a ( k/,o_l,u z )dx> ± dQdu z (3.31) 

S 0 0 0 

where 


6 iS (k / ,u 1 ,u z )3^(k / ,u 1 ,u z ) + ^(k / ,u 1 ,-u z ) 

Substituting Equation (3.26) into Equation (3.32) yields, after expanding terms and simplifying, 

f. k / 


(3.32) 


6* (k / ,o 1 ,o z ) = — -•! (t>i,o s ) + A a (o ± ,-o z )) 

~ K / l ^ 


2 —K s E' r (k/ ) + ^-( A a (oj.,o z ) - A s (o ± ,-o z )) 


Q 


(3.33) 


where 


*>- S- 

2m,. 


dfs o 
5uj_ 


(3.34) 


Using Equation (3.33) in Equation (3.31) will give a relation for the current density. However, values for 
the parameters A s (o ± ,o z ) must first be determined. By their definition of Equation (3.28), they are 
functions of the boundary values of the functions § s . 


3.3 Incorporation of Boundary Conditions 

By the definition of the problem, the plasma at the surface z = 0 is against the solid body of the 
hypersonic vehicle. Thus, reflective boundary conditions prevail from which one can write 

<l>*(0,u 1 ,u z ) = <|> a (0,u 1 ,-u z ) (3.35) 

However, at the boundary z = L , the plasma is open to the atmosphere so diffusive conditions prevail 
which give 
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<t» 5 (Z,,Oj.,-o z ) = 0 


(3.36) 


With these assignments, Equation (3.33) becomes 

(k / ,o 1 ,o z ) = n jl K 2 (Au ± ,o z )(-l) / ) 

-2^-K s E' r (k /} -^ (I,ui,i) 2 )(-l) / | 


(3.37) 


At this point, the values of the distribution functions at the boundaries, § s (0,Ox,o z ) and § s (L,o ± ,o z ) 

still remain unknown. In order to determine these values, a self-constancy condition can be applied. 
Hence, using Equations (3.35) and (3.36) with Equation (3.32) gives 

0 J (O,o 1 ,o z ) = 2f s (O,o 1 ,o z ) (3.38) 

and 

O s (^ 0 -L’°z) = ( t ) j( L ’ U -L> U z) (3.39) 

But by Equation (3.20), 

00 

O i (0,o_L,o z )= £ <M k /’ u -L’°z) (3.40) 

/= — oo 

and 

00 

0 5 (^,o_l,o z )= Yj ® 5 (k/,oj.,o 2 )(-1) / 

/ = — oo 

Hence, form these relations, one has 

00 

X < M K /’ O -L’ O 4 = 2( M 0 ’ O -L>°z) 

/=— 00 

and 

00 

£ ^s( K /’°J-»°z)( -1 ) /=< I ) 4 Z '’ 0 J-» 0 2) ( 3 - 43 ) 

/=— oo 


(3.41) 


(3.42) 


Equation (3.37) can now be used on the left sides of Equations (3.42) and (3.43) to give two equations in 
the two unknowns § s (0,o l5 o z ) and 4> 9 (Z,Ox,o z ) . To this end, substituting Equation (3.37) into 

Equation (3.42) and solving for (0,Ox,o z ) , 


<^(0,d ± ,o z ) = -£ 

/=- 00 




v Q ?- k / 2 7 


V°z J 


M*fa)+| ^ (3.44) 
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where the terms that contain the factors K; vanish, i.e., 


y iKl =o y 


(3.45) 


Similarly, using Equation (3.37) in Equation (3.43) and solving for <j) s (L,u^,u z ) yields 


<|> 5 (4oi,u z ) = -y 

/=-oo 


00 ( . 'Xl 

l 


yCl^-Kj j 




^(-1) / ^(k / ) + ^J^(I,o 1 ,o z )|| (3.46) 




At this point, it will facilitate further evaluations to simplify Equations (3.44) and (3.46) by 
analytically performing summations over the parameters / where possible. Hence, using the summations 


” ' in, ' 


z 

/=- 00 


V q 2 _k / 2 J 


iL cot (C1 s L), ^ 

/=-00 


" ' in, ' 


V q ? _k / 2 ) 


(-l/=iZ:csc(Q s Z) (3.47) 


allow Equations (3.44) and (3.46) to be written, 


<|> s ( 0 ’OL> o z) = -2 
/=-00 


iQ- 


V i \ 


v q !- k / 2 j 


V°Z J 


K s E r (k, ) - (L, oj.,0, )csc(n,Z) (3.48) 


and 




iQ , 


V , A 


/= — 00 K / J 


V°Z J 


^(- l ) / E ) ? ( K / )-%( Z , o 1 , o z ) cot ( Q , Z ) ( 3 . 49 ) 


Finally, the solution of these relations for c|) s (0,Ox,o z ) and 4 > s (Z,Oj_,u 2 ) is straightforward; after using 
some trigonometric identities, integral 


^(l,o 1 ,o z ) = (exp(2iQ,Z)-l) £ C,(/)(-l/ £*(*,) 

/ = — oo 


(3.50) 


^(0,o l5 o z )= y C,(/)^(K / )((-l) / exp(2iQ i Z)-l) 


/=-00 


( 3 . 51 ) 


where 



vn 

ko.) 


K. 


(3.52) 


Equations (3.50) and (3.51) give the sought-after expressions for the boundary values of <|> s needed in the 
evaluation of Equation (3.37). 
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3.4 Calculation of the Current Density 

Using Equation (3.2) in Equation (3.34) and evaluating Equation (3.52), one has 




V 1 

rn s l»Ts 


\ 2 




V.O* / 


z'Q* 


Q2- K/ 2 J 


/*o(N) 


(3.53) 


One can now substitute Equations (3.50) and (3.51) into Equation (3.37) and this intermediate result into 
Equation (3.31), upon remembering the definition E R (k/ ) = E R (ic/)exp(-i0) , to finally obtain the 
Fourier Transform of the current density within the plasma layer, viz, 


/ \ / ! \2 oo In oo 


I-^I J j J»l/.o(|C|)S- 2 / 


Q c 


V J A 


+ 2 i 


k/ 


u 7i 




0 0 0 


I'L l Q ?- K / 2 j 


7 


Er{*i) 


Q 2 -k 2 j 


r 


y L 'r=- oo 


/Q, 


Q 2 -k 2 ; 


exp(in, s Z,)-l| 


^i-K/ 




;Q, 


r=-»\ n s K f A Uz 7 


(- 1 ) 


(3.54) 


Using Equation (3.2) once again and performing the Uj_ and 0 integrations gives 


■M K /) = X e ^ 


( i Vi \ l/1 






. \°° / n A 

1 I f Oi 

ex P f- 

JJ t J 


| 2 1 
II 


Q c 


V j A 


v Q 2 -K 2 


+2 


K/ 


V Q 2 K 2 J 


tIz 


Q c 


V 1 ^ 


7 ! V 


Q^-k, 


L Jrt^X^s ~ K r 

i|exp(2i£J s L)-l)(-l)' £ 


/ A° z 7 

V 


Er{*i) 


\Er{ k i' )| (-1) exp(i'Q s Z)-l 
V 




Q 2 - K 2 , 


K/ i 


(3.55) 


At this point, it facilitates further calculation to change variables using the prescription 


U 7S 


which gives 


=Q J (w s ) = - 


™ S »T S 


Incorporating Equations (3.56) and (3.57) into Equation (3.55) yields 


(3.56) 


(3.57) 


Jr(k / ) = g(k / )Er(k / )- Y, a(K/,K / -)i 7 f(K / -) 

l'=-ao 


(3.58) 


where 
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a ( K l) = ^, 2ie s n s 


( 1 V 1 \ l/2 


\ m sj 


\Kj 


f 1 Yf Qj(w,) 
VTs Jj Q ?( w .s)- K / 2 


f 1 ' 


\ W *J 


exp(-w|)cfw, (3.59) 


is the single mode conductivity and 


a( K /,K/-) = ^2e 2 n, 


\ m s j 


1 

\nj 


1 

V d ts J 


I 


m 




J o l[( Q ? M~ K / 2 )( Q ? K ) - K / 2 ) 

~ j(l-(~l/ exp (**, (w, )z)j-(l-exp(2iQ iS (w* )z)) 






(3.60) 


( n *("i) _K / )( Q ?( w s)-*r)v w s 

is the multi-mode conductivity. Thus, the current density anywhere within the magnetized plasma layer 
due to the incident electromagnetic wave is given by 

00 00 00 

Jr{ z )= X CT ( K /)^( K /) ex P( z ' K / z )- X X CT ( K /’ K /')^( K /') ex P( ZK / z ) ( 3 - 61 ) 

/=— oo /=— oo/'=— 00 


At this point, contact must be made with the current density entering into Maxwell’s equations. 


3.5 Connecting the Current Density in the Plasma Layer with the Maxwell Equations 

One first needs to relate the auxiliary complex electric field E R (z) = E x (z) + iE y (z) of the incident 

wave to the governing Maxwell Equations. Returning to Equations (2.3) and (2.4) and employing 
Equations (3.5) and (3.6) gives in terms of the component fields 

~^x+^y = ik(B x x+ B y y) (3.62) 

and 

- ^ x + y = -ik(E x x+E y y) + ^-( j x x +j y y) (3.63) 

Multiplying these relations by the auxiliary vector x + iy , as used earlier, gives, respectively, 


and 


^--kB R =0 (3.64) 

oz 


8Br 

8z 


4- kEv 


.4jr . 
~ l — Jr 


c 


(3.65) 
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where (z) = 5 x (z) + iB y ( z ) and j R (z) = / x (z)+i/ > ,(z). From these two relations, the corresponding 

wave equation for E R can be derived in the usual way, 


d 2 E R 2 4n . 

— — — — + k 2 E R = -ik — j R 
OZ z c 


(3.66) 


As the equations remain amenable in the Fourier domain, the transform of Equation (3.66) is now needed 
to be used with Equation (3.58): thus, applying 


1 f 

£* ( K /) = yJ E R ( z )exp(-iKiz)dz 


(3.67) 


by multiplying Equation (3.66) by exp(-?K/z) , integrating by parts, and using Equation (3.64) yields 

r Y (k; r ( Ki ) = ( K , ) (3.68) 


Finally, substituting Equation (3.58) into Equation (3.68) and rearranging terms gives 


Kj - A (K() k («,)= 4 

c ) L 


,. (M^X-iZ-Mo)^ .,4* ^ \ 

+ZK, — — ~ lk — 2^ CT ( K /’ K /') £ «( K r) 


(3.69) 


This equation connects the fields on the surface boundaries of the plasma layer to the Fourier transform of 
the fields within the plasma. This relation can be put into a more familiar form simply by dividing 
through by k 2 using Equation (3.57) and rearranging some factors to obtain 


J- g ( K /) £*(*/)+ X s { K i’ K r)E R {Kr) 


. k ,(e r (l)(- if - e r (o) b r (l)(-i) 1 - b r (o) 


(3.70) 


s 03 o (“ + iv s ) - »Ts K f w i 


(3.71) 
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^co? 


l-(-l)^ expf i 


co-co^ +IV S 




j 2 2 

o (g> — “os + iv s ) - ^Ts K j W i l[( ra - ®cs + ™ s ) - U^K jy% 


f I !® ^CS ~^~ Ws ^ 

UftWj Jj 


(3.72) 


U7EW, (-1)' 


i+r 


® ®CS "^"^S ^ 7 j K / W 5 


1-exp 


/ III 


where the plasma frequency to /;v of the s-th charge species is given by 



4ne}n s 

m s 


(3.73) 


Equation (3.70) forms the set of equations for the transform of the electric field within the plasma layer 
that will be considered in what is to follow. 

Unfortunately, the evaluation of the integrals indicated in Equations (3.71) and (3.72) cannot be 
completed in its entirety. Although Equation (3.71) does indeed possess an analytical solution, its form 
makes impossible the analytical evaluation of the electric field from the transformed field E r {ki) in 

Equation (3.70). However, in the limits of weak spatial dispersion, the equations lend themselves to 
analytical evaluation. Of course, it must be established that the re-entry plasma can be treated as one 
possessing weak spatial dispersion. Again, it must be kept in mind that what is needed here are 
expressions for the reflection, transmission, and absorption coefficients at the z = 0 and z = L boundaries 
of the plasma layer. Thus, the transformation of the fields E R (k/) on the left side of Equation (3.70) 
must be evaluated at these surfaces. This is the goal of the following development. 


4.0 Evaluation in the Limit of Weak Spatial Dispersion 

4.1 The Limit of Weak Spatial Dispersion 

Equation (3.71) can be rewritten as 

Vrc 5 ® 0 1 “ »Ts K i ™$/ (« -©c + iv s ) 

Letting (to - co cs + iv s )~(co- co cs ) 2 + vf , consider the terms within the denominator of Equation (4. 1) in 
the case where 






«1 


which can be re-expressed as 
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(4.2) 


z((co-co c , s ) 2 +v 2 ) lWs 

Now, the integral over the variable is rapidly attenuated due to the presence of the exponential 
function. Hence, only small values of this variable need be considered. Take the value of this parameter to 
be on the order of unity. There then remains the range of values of the mode number /. In order to 
incorporate the long range effects within the plasma (for spatial dispersion effects), small values of this 
parameter need be used; thus, take / = 1 . Hence, the condition of Equation (4.2) can finally be written as 


u 


Ts 


/^((O-CD^f+vf) 


<1 


At resonance where co ~ co CiS , this condition becomes 


(4.3) 


v Ts < Lv s (4.4) 

which defines the region of weak spatial dispersion for this problem. The range over which this inequality 
holds must now be examined. 

4.2 Weak Dispersion and a Re-entry Plasma 

Since both the thermal velocity and collision frequency are functions of plasma temperature, it will be 
instructive to determine over what temperature ranges Equation (4.4) holds for nominal thicknesses and 
electron number densities. Figures 3 and 4 show the region of applicability of the constraint of Equation 
(4.4) for a one-component electron plasma of thickness L = 30 cm . As seen from these plots, the 
applicability region becomes smaller over a range of temperatures as the electron concentration decreases. 
Thus, for n e > 10 10 / cm 3 , the approximation of weak spatial dispersion can be expected to hold for 
nominal re-entry temperatures that are typically T e > 3000 K . The case in which n e < 10 10 / cm 3 requires 
one to consider the case of strong spatial dispersion. 

In the case of weak spatial dispersion for this one-component plasma {s = e). Equation (4.1) reduces 
to 


e(K/)*l 


2 <°r exp(-wj) ^ 1 < 

yfn co j (co-co ce + /v e ) e co (co-co ce + rv e ) 


(4.5) 


The analysis of Equation (3.72) proceeds along the same lines where, it must be added, that the 
exponential factors oscillate away to zero. Also, with Equation (4.4) prevailing, the terms in Equation 
(3.72) can be neglected in this approximation. In this event, Equation (3.70) can be written, using 
Equation (4.5), 

(k 3 -A: 2 + A: 2 ^(( O ))^(k / ) = /^-(^(Z)(-1) , -^(0)) + |(s r (Z)(-1) / -S r (0)) (4.6) 


where 
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(4.7) 


j/ \ 1 

A \ (ii ) = —( ■ \ 

CO (0O-CO ce + 2V e ) 

From Equation (4.6), expressions for the boundary values E R (0) and E r (L) can be obtained from 

which incident, reflected, and transmitted fields can be related finally giving rise to the associated 
coefficients. Thus, defining 

D[((o) = kJ - k 2 + k 2 A((o) (4.8) 

one has 

m°)- £ (4 - 9) 

/=- OO /=- 00 * V 10 / X ' 

and 

M*)- £ )(->)' = £ (410) 

/=- oo /=-oo / v 04 / V J 

where, as noted earlier, the term with the coefficient Kj sums to zero. 

At this point, one needs to relate the magnetic fields (0) and B R (Z) to the prevailing electric 
fields on the boundaries so as to obtain from Equations (4.9) and (4.10) two equations in the two 
unknowns E R (0) and E R (L) . From these, the reflection and transmission coefficients will be obtained. 
To this end, from Equation (3.64), 

, . 1 8E r (z) 

S *«= F -Sr (4 - n) 

As for the fields above the plasma layer at z = L , one can write 

E r (z) = Eq exp^A:(Z -z)j + exp[-ik(L -z)) (4-12) 

where Eq is the incident electric field of the wave upon the plasma layer and -ER e fl is the field reflected 
from the layer. Using Equation (4.12) in Equation (4.1 1) gives 

Br (z) = -i [£o exp(i£(T - z)) - £ Refl exp(-/A:(T - z))] (4.13) 

Hence, on the boundary z = L , 

E R (L) = E 0 + E Re fl , B R (L) = -i(E 0 -E Re ft) (4.14) 

Similarly, for the transmitted field into the region z < 0 , 

Er (z) = Ej rans cxp(z’&(Z, — z)) (4.15) 

giving at the boundary z = 0 
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Er (0) — ^jrans Cxp(?Ai), B ft ( 0 ) — //^Trans ^ P ( ) 

Substituting Equations (4.14) and (4.16) into Equations (4.9) and (4.10), defining 


(4.16) 



Figures 5 to 7 display the results of the calculation of these relations for three electron concentrations 
for a plasma of thickness L = 30 cm and T e = 3000 K [Sharkey, 2004] upon which an electromagnetic 
wave of frequency / = 2.0 GHz is incident. As Figs. 5 and 6 show, the reflection coefficient (R) remains 

at 1.0 and the transmission coefficient (7) remains at 0.0 until the magnitude of the applied magnetic field 
goes above 500 Gauss. The absorption coefficient (A) also begins to rise (clearly shown in Fig. 6). After 
this point, the coefficients rapidly oscillate and tend to their limits R — » 0.0 , T — » 1 .0 , and A — » 0.0 . For 
n e = l.Ox 10 11 /cm 3 , these limits are quickly approached. At n e = l.Ox 10 12 /cm 3 , oscillations of the 
values of R and T remain beyond the magnitude of 10,000 Gauss although the values are clearly 
separated. However, Fig. 7 for the case of n e = l.Ox 10 13 /cm 3 shows that the oscillations remain mixed 
together where the values of these coefficients exchange local maxima and minima. 


NASA/TM— 2009-21 6096 


18 



5.0 Conclusion 


It is thus demonstrated that a magnetic field of nominal magnitude applied to the plasma makes the 
plasma transparent to frequencies smaller than the plasma frequency thus substantiating the magnetic 
windows concept. Magnetic fields on the order of 4 kGauss will render a plasma transparent at 2 GHz 
with an electron concentration of 10 12 electrons/cm 3 . As the plasma density evolves throughout re-entry, 
the reflection, transmission, and absorption coefficients can oscillate over large ranges for a fixed applied 
magnetic field. Hence, a variable magnetic field controlled, e.g, by the value of received power, must be 
considered. Also, the case of strong spatial dispersion as well as the situation intermediate to weak and 
strong spatial dispersion must yet be considered. 

This example of a homogeneous plasma in a homogeneous magnetic field is certainly an idealization 
of what would exist in reality. However, these results provide impetus for further work with more realistic 
situations. Although Halbach magnets can be used to generate the applied magnetic field, some spatial 
variation of the field should be introduced. The variations inherent in the plasma thickness, temperature, 
and density must be addressed in the specification and design of a closed loop control system to adjust the 
external magnetic field to maintain transparency. Finally, account of the other components of the plasma 
must be made. In fact, with the surface ablation that occurs during re-entry, one should consider a dusty 
plasma, the kinetics of which can be vastly different from the simple one-component electron plasma 
considered here. 
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Appendix A. — A Model of the Effective Collision Frequency 

For purposes of the discussion given here, consider a two component plasma made up of electrons 
and ions in an electric field E . The Boltzmann equation including the associated collision integral for an 
electron in a two-component plasma is given by 

^+o-V r -/+-5-V c / + S coll =0 (Al.l) 

dt m 

where the complete collision integral is 

Scoll = {{?(«, e) M [/(o)F(6 1 )-/(o')F(oi)]j 3 oJQ (A1.2) 

In the case where the collision is limited to the scattering of an electron from an ion, q(u,Q) is the 
differential cross-section of electron scattering where 0 is the scattering angle, u = |o - t5j | is the relative 
velocity of the electron with respect to the ion, i.e., o and 0] are the velocities, respectively, of the 
electron and ion after collision whereas 0' and o'i are those before the collision. Additionally, is 

the velocity distribution of the ions which is usually taken to be given by the Maxwell distribution of the 
form of Equation (3.2), and dQ. = sin9d0d(j) is the differential scattering angle. 

Simplifications will now be introduced using what is known about the electrons in the plasma. First, it 
is assumed that the thermal velocity of the electrons and ions are much greater than the associated 
directed velocities in the electric field. In this isotropic plasma case (i.e., where the magnetic field 5 = 0), 
the spatial gradient of the distribution /is directed along the z axis, parallel to E . Hence, one can treat the 
“directional” part of/as a perturbation in velocity space. The distribution can thus be expanded into zero- 
order spherical polynomials, i.e., Legendre polynomials 5/ (cos 9) where 9 is the angle between E and 
o , viz., 


Writing 


/ = /(r,o,t) = ^5/(cos9)//(r,o,t) 

1=0 


^ 8f E sin 2 9 8f 

E- V 6 / = 5cos9— + — -r- 

5u o o{cosv) 


(A1.3) 


(A 1.4) 


and using Equations (A1 .3) and (A1 .4) in Equation (A1 .1) give the following system of equations for / 

S/o o dfi eE d / z \ 

-^- + — — + — (o 2 /i )+5 0 =0 

dt 3 dt 3mv 2 Su' ^ 


dt 


+ 0 


8A + 1VL 

8z 5 8z 


+ 


2 A(„3 /2 )] + s 1=0 

Su 5u 3 So' ' I 


m 


(A1.5) 


A 

8t 


+ 0 


3 8z 3 8z 


+ - 


eE 


m 


2 8 
— o — 
3 So 


'A 


+ *=° 


7o 4 So 
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where 


Si = 


2/ + 1 
471 

2/ + 1 
471 


j>,(cosft)&/© 

m«(- ,6)u/}(cOS&) 


00 00 ] 

y fl'Pj' (cosiT).F — E* y fl'Pi' (cos ft") d^u\dQ.d® 

./'= 0 /' =0 J 


= - 7*7/7) (cos 0 ))j 3 uc?Q 


(A1.6) 


In arriving at the result of Equation (A1.6), the integration over J© = sinBc/Bcftp made use of the fact 
that 


cos O' = cos 0 cos 0 ' + sin 0 sin 0 ' cos 9 

and the addition theorem for Legendre polynomials was employed. 

The chain of equations of Equation (A1 .5) can be terminated at the second one if the perturbation 
component f 2 can be neglected compared to the fundamental component /o , i.e., 


5/o 1 8 

— »— — 
So o 3 So 


(o 3 / 2 ) 


(A 1.7) 


Before Equation (A1 .7) can be established, one now needs to consider the expressions for the scattering 
perturbations Sq, S\, and S 2 . To this end, considerations will only be made of elastic collisions; in this 
case, Equation (A1 . 6 ) gives for S\ 


,0)m(/i(o)f(o,)-^(cos0)/i(o')F(o'i))^ 3 o 1 JQ (A1.8) 


Assuming that the energy of the electron changes only slightly after collision with the ion, one has that 
|o'| « |o| and |o'i | « | 6 ] | where |o| » |oi | . In this instance, Equation (A1 . 8 ) yields 

Si =/i(o)JJ?(m,0)oF(oi)(1-cos0)c1 3 o 1 ^Q =v 1 (o)/ 1 (o) (A1.9) 

where the velocity dependent collision frequency is defined by 

vi(«) = Arj on oj< 7 (o, 0 )(l-cos 0 )rfQ, AT ion = jV(oi)*/ 3 Oi (A1.10) 

Equation (A 1.9) begins to form the basis of the Krook model introduced in Equation (2.7) where a 
velocity independent collision frequency v 5 was used. It is the central purpose of this Appendix to obtain 
an expression of this velocity independent collision frequency from collision theory. Before this can be 
accomplished, however, it remains to establish the prevailing conditions that allow the termination of the 
system of equations of Equation (A1.5). 

One can similarly show that 
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S 2 = JJ^(«,0)M(/2(6)F(6i)-^(cos0)/ 2 (o')F(oi))rf 3 o 1 £/Q 
= /2(O)||?(m,0)oF(o 1 ) 

= V 2 (u)/ 2 (u) 


^ 3cos20 + l 




(Al.ll) 


where 


v 2 (M) = A io no|?(u,0) 


1- 


3cos20 + 1 


\dCl, N ion =Jf(ui)^ 3 ui 


(A1.12) 


Thus, the quantities v\ (o) and v 2 (o) are of the same order of magnitude and, for purposes of the 
discussion to follow, v(o) = v\ (o) « v 2 (o) . Finally, an expression for Sq needs to be secured. The 
calculation is trivial upon using Equation (A1 .6) and noting the fact that energy exchange during 
collisions does not occur with the approximations used above; hence, to first order, Sq = 0 . Finally, 
specializing to the special case of spatial homogeneity, S/Sz = 0 in Equations (A1.5). Given the 
developments above, these relations reduce to 


8t 


eE 8 
3 mo 2 So 



= 0 


(A1.13) 


8f\ | eE a/ 0 | 
8t m So 


v(o)yi =0 


Qh + ^E_ 
8t m 


2 8 
— o — 
3 So 


( i A 




fl + 


3 8 

7o 4 So 


(u 4 / 3 ) 



(A1.14) 

(A1.15) 


Consider now the steady state situation where /J,/ 2 ~exp(-zcot) and hence 8f/8t — zco/j and 
S/ 2 /St — iobfi . Using these in Equations (A1.14) and(A1.15), dropping the second term within the 
parentheses of Equation (A 1.1 5) and solving for |/ 2 | gives 


e 2 F 2 


m 1 - 


! (co 2 + v 2 ) 

Hence, the condition of Equation (A 1.7) becomes 

/o» 


o 


So 


1 V6 

o So 


(A1.16) 


e 2 E 2 

f 1 ^ 

8 | 


(ft) 2 + V 2 ) 

V u J 

dv 

l 30 ; 


If Equation (A 1.1 7) holds, the original system of equations of Equation (A 1.5) reduce to two: 

i + £4(^H 


(A1.17) 


(A1.18) 
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Finally, if the electric field is weak enough where the second term of Equation (A1.18) can be taken as a 
first order perturbation, the symmetric portion of the distribution / retains its Maxwellian nature and the 
solution to the problem reduces to (remembering Equation (A1 .3)) 


/ = /o + /icos9 


d/i ( eE df 0 | 
8t m e So 


v(o)/i =0 


(A 1.20) 


In what is to follow, it is found advantageous to give the perturbation f a vector character by writing 
Equation (A 1.20) as 


/-/o + 


7i-g 

o 


(A1.21) 


g/i | eE a/p | 
dt m e do 


v(u)/i=0 


(A 1.22) 


where f\ is given the direction of E . 

The total current induced in the plasma by the free electrons is given by 

j =en e \of(v)d 3 v (A1.23) 

Using Equation (A1 .21) and converting the u -integration into one in spherical coordinates gives, noting 
the isotropy of fo , 


oo 271 

j = 2nen e j J o 3 /i sin0c?0^o 
0 0 


Now, from Equation (A1 .22), one has in the steady state 


eE d/o 
m e do 
v(o)-fco 


Using Equation (3.2) in this result and substituting into Equation (A1 .24) yields 


-r %en a E 
J = 




7rm„ 


l 


v(x)x 4 exp^-x 2 ) “Sx^exp^-x 2 ) 


o 


co 2 + v\x) 


dx + ico J 


co 2 +v 2 (x) 


dx 


(A 1.24) 


(A1.25) 


(A 1.26) 


where x = o^m e /2A:g7’ e . In the event that the collision frequency is independent of the velocity, i.e., 
v( x) = v = const. , the integrals in Equation (A1 .26) can be performed and give the result 
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V 


1 


(A 1.27) 



v ft) 2 + v 2 


+ i ft) 


ft) 2 + v 2 J 


This result is identical to the one that is obtained by using elementary considerations. That is, using the 
concept of an effective collision frequency v eff , one simply uses Newton’s Law of motion for an electron 
in an electric field and writes 


m„ 


d u 
dt 


= eE- w e v e ff o 


(A1.28) 


Here, the effective collision frequency enters as a friction term. Again, considering the steady state and 
solving for the velocity, one gets for the associated current 


J = en e o = 


e 2 n e E 

m„ 


v eff 
CO 2 + v[ 


- + ZOO 


eff 


C0 2 +V 2 ff J 


(A 1.29) 


In order to reconcile these two approaches, one considers the limiting case where co 2 » v 2 ff ,v 2 (x) . 
Taking Equations (A1 .26) and (A1 .29) in this limit, equating the two results and solving for v e ff gives in 
terms of the velocity dependent collision frequency 


g ou 

v e ff =^-y=|v(x)x 4 exp^-x 2 )cfr 


(A 1.30) 


It is now necessary to find an expression for v(x) defined by Equation (A1.10), given the scattering 

situation assumed here. That is, elastic scattering of a fast moving electron from an essentially stationary 
heavy ion. In this Coulomb scattering case, one must use the Rutherford scattering formula for the cross- 
section required in Equation (A1.10), viz., 


9(0, 0) = 


' e 2 * 


v 2w e u 2 j 


sin 


4(0, 


(A1.31) 


Substituting this expression into Equation (A1 . 1 0) and performing the azimuthal integration, 


( e 2 If 

v(u)=2™ i “,u|l— j j J 


1-COS0 . ... 
-sin0J0 


27t«i 0 


( e 4 > 

v w 2 o 3 J 


sin 


4/0 


(A1.32) 


In 1 + cot 2 




J) 


where ©min is the minimum scattering angle which is related to the maximum value of the impact 
parameter b max , 


tan 


e 


2 


BJgU 2 6 max 


(A1.33) 
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Here, the impact parameter is determined by the fact that in a plasma where the interaction between the 
ions and electrons is through a Coulomb field only, distances on the order of the Debye radius r D cannot 
be exceeded since this is the maximum distance at which substantial interaction occurs between the 
electron and ion; at distances greater than this, the field of the ion decreases exponentially. Hence, one 
can write 


^max ~ r D ~ 


ion^e 


vl/2 


47ie 2 «i on (J^bI ion + ^B^e) 


k B T e 
Sne 2 ) % 


\l/2 


ion J 


(A 1.34) 


where the last result issues from the assumption that 7j on « T e . For a typical re-entry plasma, 

r D ~ 10 -4 cm . Given the relative values of the parameters involved, one can write from Equation (A1.33), 


0 min — 2 tan 


m e° r D 


2e 2 

m e v 2 r D 


(A1.35) 


Substituting this result into Equation (A 1.32) and series expanding the cot 2 function finally gives 


v(u) = 27I«i 0 


f .4 \ 


2ii3 


V W|U j 


In 


t | r 2 m 2 u 4 ^ 


(A 1.36) 


Finally, applying the substitution o = x^2k B T e /m e in Equation (A1 .36) and using this result in Equation 
(A 1.30) yields 


v eff = 


3^ 


(27t«ion) 




\3/2 oo 


m„ 


\ m e A^k B T e 


J jcln(l + ,4 2 jc 4 jexp^-jc 2 'jdx 


3yfn 


2nni 0n e 4 


m,. 


\3/2 


mi 


2k B T e j 


1 . 
—sin 


f i A 




7i -2 si 


\Ajj 


■cos 


f \ vmi 


(A1.37) 


ci 






where A = 2k B T e r D l e 2 and ci(---) and si(---) are the cosine and sine integrals, respectively. For a 

typical re-entry plasma, A ~ 10 2 » 1 . Thus, expanding the functions within the brackets of Equation 
(A 1.37) in an ascending series, the sine terms are negligible and one is left with 

cos(l/^)ci(l/^) « -In A + y where y is Euler’s constant, y « 0.577 . Hence, Equation (A1.37) reduces to 


v e ff : 


3yfn 


f r. JV ^ /2 f 

In 


27t« ion e 4 Y m e 


mi 


2k B T e j 


\.\2k B T e r D 


(A 1.3 8) 


This is the expression for the effective collision frequency used in this work. To make contact with the 
notation used in the text, one has v v = v e ff where s = e . It must be remembered that this formulation 
only accounts for elastic collisions between electrons and heavy ions. Of course, other scattering 
processes can occur such as elastic and inelastic collisions with molecules, collisions with dust grains in 
which charge transfer can also be attendant with impact, etc. 
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B 


Incident EM Field 
E = E x x + E y y 

B=B x x+B y y 


Figure 1. — The flowing plasma and its various fields and velocities. 
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Figure 2. — Electromagnetic and kinetic equations for a plasma in an external 
magnetic field and an incident electromagnetic wave. 
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Figure 3. — Region of applicability for weak spatial dispersion approximation for n e =10 10 electrons/cm 3 and L=30 cm. 
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Figure 5. — Reflection, transmission, and absorption coefficients versus applied magnetic field strength for 
fie^O 11 electrons/cm 3 . 
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Figure 6. — Reflection, transmission, and absorption coefficients versus applied magnetic field strength for 
r? e =10 12 electrons/cm 3 . 
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Figure 7. — Reflection, Transmission, and Absorption Coefficients Versus Applied Magnetic Field Strength for n e = 
electrons/cm 3 . 
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